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1 . Introduction 

Heymsfield et al. (2001) described the structure of 
Hurricane Bonnie on August 23, 1998, using multiple 
observations from the NASA Convection And Moisture 
Experiment (CAMEX-3). On that day, Bonnie's structure 
wasi highly asymmetric with strong low-level reflectivities 
on the eastern side of the storm and little precipitation on 
the western side. Observations from the NASA ER-2 
Doppler radar indicated isolated deep convective towers 
on the eastern side of the storm with updrafts greater 
than 10 m s \ Heymsfield et al. (2001) developed a 
conceptual model for the evolution of the convective 
towers in which the updrafts formed near the top of the 
boundary layer on the southern side of the eyewall and 
grew progressively taller while moving around to the 
northern side. In other words, the convection was viewed 
in terms of rising plumes of air that were initiated to the 
south and reached the upper troposphere to the north of 
the eye. 

this study examines a high-resoiution simulation of 
Hurricane Bonnie. Results from the simulation will be 
compared to the conceptual model of Heymsfield et al. 
(2001) to determine the extent to which this conceptual 
model explains vertical motions and precipitation grovirth 
in the eyewall. 

2. Methodology 

The model used in this study is the PSU-NCAR non- 
hydrostatic fifth generation mesoscale model (MM5 
V3.4). The simulation uses four levels of grid nesting 
starting at 36 km and then nesting to 12, 6, and 2 km. 
The simulation is started at 1200 UTC 22 August 1998 
and run for 36 hours with only the 36- and 12-km grids 
active, with model output saved every hour. A one-way 
nest is used to conduct a higher resolution simulation on 
the 6- and 2-km grids starting at 6 h into the simuiation. 
The higher resoiution grids are run untii 36 h with output 
every 15 min. A 6-h period from 24-30 h is run with 
output every 3 min. The 2 km grid was moved hourly to 
keep it centered on the storm. Physics options inciude a 
modified version of the Blackadar planetary boundary 
layer scheme in which parameterization of the surface 
roughness parameters follows Garrett (1998) and 
Psigowski and Moore (2001). Cloud processes are 
represented by the Grell cumulus parameterization 
scheme (on the 36- and 12-km grids only) and the 
Goddard cloud microphysics. Radiative processes are 
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calculated every five minutes. 

Initial and boundary conditions were obtained from 
12-hourly global ECMWF analyses archived at NCAR. 
Because the large-scale analysis does not contain an 
adequate representation of the initial hurricane vortex, a 
bogusing technique using four-dimensionai variationai 
data assimilation is used (Pu and Braun 2001). The 
assimilation was performed on the 36-km grid only. 

The storm center is determined, as in Braun (2002), 
at every model output time using the pressure field at the 
lowest model. The horizontal distribution of pressure is 
used to determine an approximate geometric center, or 
centroid, of the pressure field. The location of the 
minimum pressure is used as a first guess for the center. 
A variational approach is then used that adjusts the 
location of the center until the azimuthal variance of the 
pressure fieid at all radii between the center and the 
outer portion of the eyewall (~ 65 km) is minimized. This 
methodology works well not only at identifying the 
centroid of the pressure field but also the approximate 
centroid of the ring of strong tangential winds and 
vorticity. Storm motion is then computed from the 
identified center locations. 

3. Wavenumber 1 asymmetry in the eyewall 

In this section, we investigate the time-averaged 
vertical motion and precipitation fields and the 
mechanisms that force the wavenumber 1 asymmetry. 
Figure 1 shows the 6-h averaged fields of vertical motion 
and total precipitation mixing ratio (sum of rain, snow, 
and graupel) at four levels from the top of the boundary 
layer to the upper troposphere. The precipitation 
contours show that the maximum precipitation occurs on 
the northeastern side of the storm at all levels. In 
contrast, the vertical motions show some variation with 
height. At 1 km, the maximum upward motion is on the 
northeastern side of the storm just inside the region of 
maximum precipitation. At mid levels (5 and 8 km), the 
maximum upward motions occur on the southeastern 
side of the storm. At 12 km, the upward motions are 
much weaker with no clearly defined maximum. 
Downward motion is strongest at mid-to-upper levels 
and is concentrated in two areas: on the western side of 
the storm within or just outside the western portion of the 
eyewall and within the eye just inside the heavy 
precipitation associated with the eyewall. 

Many studies have examined the impact of vertical 
wind shear on the development of wavenumber 1 
asymmetries in the eyewall and have described several 
mechanisms by which this occurs. These include: (1) the 
tilting of adiabatic vortices by the shear and the 
subsequent development of upward motion in the 
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Figure 1. Time averaged fields of vertical motion (colors, intervals of 0.25 m s'’ for w<0. 0.5 m s'’ for w>0) and total 
precipitation mixing ratio (sum of rain, snow, and graupel, solid contours at intervals of 0.5 g kg'’). Fields are 6-h 
averages using output every 3 min from hours 24-30 of the simulation. 


downtilt direction as required for the flow to remain 
balanced (Jones 1995; Wang and Holland 1996); (2) the 
development of upward motion 90 degrees to the right of 
the tilt direction as a result of the interaction of the vortex 
flow with the temperature asymmetries generated by (1 ); 
and, (3) the effects of flow relative to the storm and an 
assumed balance between horizontal vorticity advection 
and vorticity stretching or compression that generally 
produces upward motion in the downshear direction 
(Bender 1997; Frank and Ritchie 2001). 

To examine mechanism (1 ), the center position from 
the time-averaged fields was determined at each height 
and the results (not shown) indicate a south-eastward tilt 
of 6-8 km between heights of 1-12 km. Figure 1 
therefore shows that the maximum upward motion at 
mid-to-upper levels occurs in the downtilt direction. This 


result suggests that the tilting of the vortex may play a 
role in the development of the asymmetry and ftjrther 
suggests that mechanism (2) does not. Mechanism (3) 
was examined by overlaying contours of horizontal 
divergence with contours of vorticity as well as wind 
vectors associated with the asymmetric wind 
(determined by subtracting out the azimuthally averaged 
tangential and radial winds). According to the vorticity 
balance argument, areas of negative (positive) vorticity 
advection should be associated with areas of 
convergence (divergence). The results (Fig. 2) indicate 
that regions with asymmetric inflow (outflow) in the 
eyewall are associated with convergence (divergence). 
The relative flow in the boundary layer (Fig. 2a) is from 
the northeast so that the maximum boundary layer 
convergence occurs on the northeastern side of the 
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eyewall, thereby producing the upward motion maximum 
there in Fig. 1a. Above the boundary layer, the relative 
flow produces low-level convergence (Fig. 2b) and 
upper-level divergence (Fig. 2c) on the southeastern 
side of the storm, thereby producing the maximum 
upward motion on that side of the eyewall at 5 and 8 km 
(Figs. 1b, 1c). Consequently, the pattern of upward 
motion at most levels can be explained by mechanism 
(3). 

4. Updraft structure and eyewall mesovortices 

Although the time-averaged vertical motion fields 
show a relatively smoothly varying asymmetric vertical 
motion pattern, the upward motion in the eyewall at any 
instant in time is typically comprised of one or more 
convective updrafts that occupy only a relatively small 
percentage of the eyewall area (as in Braun 2002). 
These updrafts are observed to form on the southern 
side of the eyewall and dissipate on the northern side, 
similar to the conceptual model of Heymsfield et al. 
(2001 ). However, in contrast to the conceptual model. 
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Figure 2. Horizontal cross sections of divergence 
(colors), vorticity (contours), and asymmetric winds at 
(a) 245 m, (b) 1.1 km, and (c) 8.2 km. The contour 
interval for divergence is 0.5x10'^ s \ with the zero 
contour highlighted by the dotted contours. For 
vorticity, the contour interval is 1x10"^ s'^ starting at 
1x10"^ s'\ The wind vector above the upper level 
comer corresponds to 5 m s'\ 

the updrafts are not starting off shallow on the southern 
side and growing to maximum height on the northern 
side, but are instead seen to extend through the 
troposphere virtually from their inception. Also, the 
updrafts are seen to first develop at mid-upper levels (3- 
10 km) on the south-southwestern side of the storm and 
then to extend downward to near the top of the boundary 
layer on the southeastern side of the storm. This pattern 
of development can be seen in Fig. 1 by the shift in the 
location of upward motion between Fig. la and Figs. 1b 
and 1c. 

To explain the mechanisms that control the timing 
and location of the updrafts, consider the distributions of 
radial velocity, vorticity, and upward motion at 1.1 km 
shown in Fig. 3. The relative flow at this level is from the 
southeast (Fig. 2b) such that the strongest inflow occurs 
on the southeastern side of the storm and outflow occurs 
on the northwestern side. In Fig. 3a, a ring of high 
vorticity is seen associated with the eyewall and within 
this ring are areas of highly concentrated vorticity, or 
mesovortices. Considering the mesovortex on the 
southeastern side, its circulation should produce relative 
outflow (inflow) on its trailing (leading) side, where the 
terms “trailing" and “leading” are defined in terms of the 
movement of the mesovortex cyclonically around the 
eyewali. The radial velocities show a region of weak 
outflow trailing the mesovortex and an area of inflow 
extending into the eye ahead of the mesovortex, 
consistent with this interpretation. Where the outflow 
meets with the general area of inflow, there is enhanced 
convergence and, in the area of vortex-induced inflow, 
there is reduced convergence. As a result, the enhanced 
convergence trailing the mesovortex leads to the 
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Figure 3. Horizontal cross sections of radial velocity at 1 .1 km and (a) vorticity and (b) vertical motion. The implied 
circulation associated with the southeastern mesovortex is indicated. The contour interval for radial velocity is 5 m s'\ 
with the zero contour highlighted by dotted contours. For vorticity, the contour interval is 1x10"^ s'^ starting at 3x10'^ s’’ 
while for vertical velocity the contour interval is 1 m s'^ starting at 1 .5 m s \ 


formation of the updraft there (Fig. 3b). 

To show that this pattern exists at other times, Fig. 4 
shows a Hovmoller diagram of the radially averaged 
absolute vertical vorticity and vertical velocity at 1.1 km. 
The vorticity field shows positive vorticity anomalies that 
rotate around the storm at least once and, in some 
cases, several times. When these positive vorticity 
anomalies move into the eastern half of the storm, 
updrafts form and slightly trail the vorticity anomalies in 
time in a given direction. Consider vorticity anomaly A in 
Fig. 4. It starts on the western side of the storm and 
initiates an updraft when it moves into the eastern side 
the storm. It subsequently rotates around the storm, 
initiates a new burst of convection on the eastern side, 
and then rotates around again to eventually form a third 
round of convection. Vorticity anomaly B is associated 
with two episodes of convection before dissipating on 
the southwestern side of the storm. In some cases, such 
as vorticity anomaly C (formed by the merger of two 
vorticity anomalies), the vorticity appears to be 
generated by the convection and then subsequently 
moves around the storm to initiate a subsequent episode 
of convection. These results suggest a strong linkage 
between eyewall mesovortices and convective updrafts. 
Updrafts form when the mesovortices move into the 
region where the relative flow supports the upward 
component of the wavenumber 1 asymmetry in vertical 
motions. 

Given this understanding that the updrafts form 
where the mesovortex interacts with the wavenumber 1 
forcing of vertical motion, we are now in a position to 
explain why the updrafts form at mid-to-upper levels 
prior to forming near the top of the boundary layer. 
Within the boundary layer, the relative asymmetric flow 


produces maximum convergence on the northeastern 
side of the eyewall (Fig. 2a) so that the interaction of the 
mesovortex with this convergence pattern will produce 
updrafts that form on the southeastern side of the storm 
and dissipate on the northwestern side. In contrast, 
above the boundary layer, the asymmetric flow produces 
low-level convergence (Fig. 2b) and upper-level 
divergence (Fig. 2c) on the southeastern side of the 
storm, allowing mesovortices to initiate updrafts on the 
south-southwestern side and later dissipate on the 
northeastern side of the storm. 

The mesovortices form as a result of the breakdown 
of the ring of vorticity within the eyewall. As shown by 
Schubert et al. (1999), a ring of vorticity is barotropically 
unstable so that small perturbations will grow rapidly 
allowing the axisymmetrization of vorticity to gradually 
change the ring of vorticity to a vortex monopole. In a 
storm with active convection, the axisymmetrization 
process does not reach this final state of a vortex 
monopole because the convection maintains the ring of 
vorticity. Consequently, mesovortices are continually 
generated as the ring of vorticity breaks down as a result 
of barotropic instability but is maintained by convection. 

5. Summary 

The numerical simulation of Bonnie suggests that 
updrafts in the eyewall are the result of the interplay 
between wavenumber 1 forcing of vertical motion and 
eyewall mesovortices. The wavenumber 1 forcing of 
vertical motion is consistent with the idea that horizontal 
vorticity advection is roughly balanced by the stretching 
or compression of vorticity such that convergence occurs 
in association with inflow into the eyewall and 
divergence occurs with outflow from the eyewall. The 
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Figure 4. Hovmoller diagram of radially averaged (20-54 km) absolute vertical vorticity (colors, contour interval of 
0.25x10'^ s'’ starting at 1.25) and vertical velocity (contour interval of 0.1 m s'’ starting at 0.6) at 1.1 km elevation. 
The plot has been extended to three revolutions around the center. Solid lines follow individual vorticity disturbances 
or mesovortices. 


mesovortices, which rotate cyclonically around the 
eyewall, contribute to upward motions by producing 
enhanced (weakened) convergence on their trailing 
(leading) sides as they move into the favorable portion of 
ttie wavenumber 1 asymmetry. 
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